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Abstract 

We present a calculation of 7r° photoproduction on the proton in manifestly Lorentz-invariant 
baryon chiral perturbation theory up to and including chiral order q 4 . With the results we analyze 
the latest ir° photoproduction data in the threshold region obtained at the Mainz Microtron. In 
the calculation of observables and the fit of the low-energy constants, we take S, P, and D waves 
into account. We compare the results for the multipoles with the corresponding single-energy 
analysis. Furthermore, we also fit the 0(q 4 ) heavy-baryon chiral perturbation theory calculation 
and compare both results. We provide predictions for several polarization observables for future 
experiments. Finally, we discuss the (3 parameter of the unitarity cusp which is related to the 
breaking of isospin symmetry. 

PACS numbers: 12.39.Fe, 13.60.Le 
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I. INTRODUCTION 



Quantum chromodynamics (QCD), the gauge theory of the strong interactions, possesses 
an approximate global symmetry due to the small masses of the u and d quarks. Exploiting 
this so-called SU(2) L x SU(2) R chiral symmetry and its spontaneous breaking to SU(2) V , 
model-independent predictions can be made for physical observables. Such a prediction 
is called a low-energy theorem (LET) and, in the case of pion photoproduction, the first 
derivation of a LET for the charged channels by Kroll and Ruderman [lj was based on elec- 
tromagnetic current conservation. For neutral pion photoproduction, the first sufficiently 
precise experimental data @, [sj] for the S'-wave electric dipole amplitude E 0+ [4| showed a se- 
rious disagreement with the prediction based on current algebra and the partially-conserved 
axial- vector current hypothesis jB], 0|. In Refs. 0, [sj, the discrepancy between experiment 
and theory was addressed with the aid of chiral perturbation theory (ChPT) which is an 
effective field theory of QCD at low energies based on chiral symmetry (see, e.g., Refs. 

I3| for an introduction). In particular, it was shown that E 0+ gets modified by certain 
nonanalytic loop contributions. Using heavy-baryon ChPT (HBChPT) (ill. Bernard et 
al. analyzed neutral pion photoproduction in Refs. (l3l-[l5l|. One motivation for introducing 
HBChPT was the fact that manifestly Lorentz-invariant (or relativistic) ChPT (RChPT) 
seemingly had a problem concerning power counting when loops containing internal nucleon 
lines come into play. In this diagram apparently has contributions which are of 

lower order than determined by the power countin g. B y choosing appropriate renormaliza- 



tion conditions, this problem was solved in Refs. [16|, |17|. The crucial difference between 
HBChPT and RChPT is that, at a given order, the latter also includes an infinite number 
of higher-order corrections. These corrections can be important as, e.g., in the case of the 
scalar nucleon form factor where one even gets the wrong analytic behavior in HBChPT 



16|| . Another example is the Fubini-Furlan-Rossetti sum rule [18j . The expansion around 



the nucleon mass shifts the pole positions in the s and u channels away from the physical pole 



positions [19|. A further, non-perturbative approach for studying pion production beyond 



the threshold region in a covariant way was developed in Ref. [20|. The method makes use of 
the chiral effective Lagrangian up to and including 0(q 3 ) and is based on the implementation 
of causality, coupled-channel unitarity, and electromagnetic gauge invariance. 

In pion photoproduction, some of the higher-order corrections turn out to be large. While 
in Ref. [l3| an 0(q 3 ) calculation obtained LETs for the P-wave amplitudes Pi and P2, the 
calculation at 0(q 4 ) [l5| gave large corrections to these LETs. As one can see from the 
numerical values of the low-energy constants (LECs), there are still important higher-order 
contributions missing. Here, we present a full one-loop 0(g 4 ) RChPT calculation. We 
find that some of these higher-order contributions are included in the relativistic case. In 
addition, we also analyze D waves and show that at order 0(g 4 ) there is another LEC, which 
mainly affects the E 2 - multipole and, through mixing, also E Q+ . We will focus on the latest 
data obtained at the Mainz Microtron (MAMI) [2l[ which give very precise results for the 
differential cross section and the photon asymmetry £ in the threshold region and which 
are, therefore, well-suited to pin down the LECs and, with that, the multipoles. 

This work is organized as follows. In Sec. II we present the formalism for neutral pion 
photoproduction. Section III gives a short introduction into the framework of ChPT. In Sec. 
IV we present and discuss our results. We give predictions for several polarization observables 



and compare them with the predictions of the Dubna-Mainz- Taipei (DMT) model [22J. We 



also show a fit of the HBChPT results to the new data to have a better comparison with 
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RChPT. Finally, we analyze the so-called f3 parameter of the unitarity cusp [23l. |24|. Section 
V closes with a short summary. Some technical details of our calculation can be found in 
the appendices. 



II. MATRIX ELEMENT 




FIG. 1: Kinematics of pion photoproduction. The solid line represents the nucleon with incoming 
momentum pi and outgoing momentum pf, respectively. The wiggly line represents the photon 
with momentum k and the dashed line denotes the pion with momentum q. The blob symbolically 
stands for all contributions to the process. 

Pion photoproduction on the nucleon is the creation of a pion from a nucleon via the 
absorption of a photon, 

N + j-> N* + ic, (1) 

where N (JV') denotes the nucleon in the initial (final) state, 7 represents the photon, and 7r 
symbolically stands for the appropriate pion. Here, we only discuss the production of a tt° 
on the proton. The kinematics of the process is depicted in Fig. [IJ The usual Mandelstam 
variables s, t, and u are defined as 

s = ( Pi + kf = { Pf + q) 2 , u = fa - q) 2 = { Pf - k) 2 , t = { Pi - p f ) 2 = {q - k) 2 , (2) 

and fulfill 

s + t + u = 2m 2 N + M 2 , (3) 

where and M n denote the nucleon mass and the pion mass, respectively. In the center- 
of-mass (cm) frame, the energies of the photon, E^, and the pion, E n , are given by 

_ W 2 -m 2 N _ W 2 + M 2 -m 2 N 
7 ~ 2W ' n " 2W ' 1 ] 

where W = <Js is the cm total energy. In the lab frame, the photon energy E l ^ h is given by 

Piab _ W 2 -m 2 N 

E -> " 2m N ■ (5) 
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The cm scattering angle G,,- between the pion three-momentum and the z-axis, defined by 
the incoming photon, can be related to the Mandelstam variable t via 

t = M*-2(E J E 7T - \k\\q\ cos (6) 

For linearly polarized photons, the differential cross section ^ in the cm frame can be 
written as 

^(e ff ,0) = 0b(e T )[i-£(e*)cos20, (7) 

where ao(Q w ) is the unpolarized cross section and ^(Q^ denotes the photon asymmetry. 
The azimuthal angle <ft is defined as the angle between the polarization vector of the photon 
and the reaction plane spanned by the nucleon and pion three-momenta. 

Using the so-called Ball amplitudes [26| , the matrix element of pion photoproduction can 
be parametrized in a Lorentz-covariant way, 

-iee,(N'7r\J»(0)\N) = e„u(p f ) fe^j ^ 

In Eq. ([8]), e M denotes the polarization vector of the photon, J M is the electromagnetic current 
operator in units of the elementary charge e > 0, and u(pi) and u(pf) are the Dirac spinors 
of the nucleon in the initial and final states, respectively. In the following, our convention 
differs slightly from Ball's original definition. We use 

Vf = 7 M 7s, V£ = 75^", 

Vi 1 = 7 5<A V£ = 75*", 

Vg = 7^,7^75, V£ = K-f-ysP", 

vf = Ki v im\ vi = k v Yi^, (9) 

with P = (pi +pf)/2. Only six of the amplitudes survive in pion photoproduction, as 
e ■ k = 0. The remaining amplitudes are related because of current conservation, 

B 1 + B 6 k-P + B 7 k-q = 0, B 2 k ■ P + B 3 k ■ q = 0, (10) 

so one ends up with only four independent structures. 

In the threshold region, pion photoproduction is commonly analyzed in terms of a multi- 
pole decomposition. The usual definition of the matrix element in the cm frame [ij is related 
to Eq. (JED via 

<*«(p/) (E By n u ^ = —XfFxi' ( n ) 
\<=i / mN 

where x% an d Xf denote initial and final Pauli spinors. In the Coulomb gauge (e = 0, 
k ■ e = 0), J 7 may be written as follows, 

JF = ia • tT\ + a ■ q 3 ■ (k x e) JF? + i<J ■ k q ■ eJ-3 + ia ■ q q ■ eT^, (12) 

where q and k denote unit vectors in the direction of the pion and the photon, respec- 
tively. Introducing x = cos = q ■ k, the so-called Chew-Goldberger-Low-Nambu (CGLN) 
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amplitudes Ti can be expressed in terms of energy-dependent multipole amplitudes as 



oo 

T x =J2 { [ lM i+ + Ei + ]Pl +1 {x) + [(I + 1)M,_ + Et_]PU(x)}, 
1=0 

oo 

oo 

1=1 

oo 

= Yl { Ml + - El + - Ml - - E i-} p !'( x )- 



(13) 



In Eq. (fl3|) . P;(x) is a Legendre polynomial of degree /, P( = dPi/dx and so on, with 
/ denoting the orbital angular momentum of the pion-nucleon system in the final state. 
The multipoles Ei± and Mi± refer to transitions caused by electric and magnetic radiation, 
respectively, and the subscript 1+ denotes the total angular momentum j = I ± 1/2 in 
the final state. These multipoles have several advantages. First of all, one can match 
their quantum numbers to those of nucleon resonances in order to analyze the excitation 
spectrum. Furthermore, while the amplitudes Ti depend on the cm angle of the reaction in 
a complicated manner, this dependence can be completely projected out in the case of the 
multipoles, as the Legendre polynomials form an orthogonal basis. The results read 27 



dx 



PiJ 7 1 — Pl+lF 2 



i 2(/ + l) 

J—riPi-x - P l+l )JF 3 + ^\(Pi ~ Pl+2^4 



Ei. 



21 + 1 
1 dx 
1 ~2A 
l + l 



21 + 3 



PlJ~l — P\-\J : 2 



Ml 



21 + 1 

1 dx 



(Pi-! ~ Pl + l)T 3 + 



I 

21 - 1 



(Pi - Pl-2)T± 



,2(1 
1 dx 
t 2[ 



P\F 1 _ Pl+lT 2 



P X T! + P X - X T % 



21 + 1 



(Pi 



i-i 



Pi+i)^ 



1 



21 + 1 



(14) 



In addition, in the threshold region one needs only few multipoles to describe physical 
observables. Traditionally, only S and P waves were used, resulting in the following simple 
expressions for the unpolarized differential cross section ctq and the photon asymmetry S 



28 



0"o(©tt) 

S(0.) 



in 

lifel 



(A + B cos Q n + C cos 2 9^) 



\q\ sin 2 
2|fc|ao(6^) 



(15) 
(16) 
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where 

A=\E 0+ \ 2 + ±(\P 2 \ 2 +\P 3 \ 2 ), 
B = 2Re(P 0+ P*), 

C=\P 1 \ 2 -±(\P 2 \ 2 + \P 3 \ 2 ). (17) 

The above-mentioned linear combinations of P waves, Pj, are defined as 

Pi = 3Pi+ + M 1+ -M!_, 
P 2 = 3P 1+ - M 1+ + Ml, 

P 3 = 2M 1+ + M!_. (18) 

Furthermore, we show predictions for polarization observables, namely, target asymmetry 
T, recoil polarization P, and beam-target asymmetries E, F, G, and H [28], 

r(©ir) = I Sm ^ {lm[P * + (^ - Ps)] + cos6Jm[P*(P 2 - P 3 )]} , 

p(e ^ } = _ | sm e {Im[K+( p 2 + p 3)] + cos ej m[ P*(P 2 + p 3 )]} , 

{|P 0+ | 2 + Re(P 3 P 2 *) + 2cos9^Re(P 1 P* + ) 

|fc|cro(©7r) 

+ cos 2 6^[|P 1 | 2 -Re(P 3 P 2 *)]}, 
ne w ) = ,'fj ^ ®; {Re[Po* + (^ - Pa)] + cose.RelP^P, - P 3 )]} , 

P(6.) = | Sm ®; {lm[P* + (P 2 + P 8 )] + cos6Jm[P*(P 2 + P 3 )]} - (19) 



In Refs. |29l. |30| . the importance of D waves was pointed out, especially of the P 2 _ multipole. 
Even though their numerical values are small, they strongly affect the extraction of other 
multipoles through interference with large P waves. The relevant formulas for the observ- 
ables are rather lengthy when the D waves are included, so we refer to Ref. [30[ for further 
details. Usually, in the threshold region the D waves are assumed to be given entirely by 
Born contributions. In the next section, a low-energy constant is discussed which explicitly 
influences P 2 _ and, because this multipole strongly mixes with E 0+ , can change the deter- 
mination of E 0+ significantly. In the threshold region, the multipoles M.i± {M. = E, M) 
are proportional to \q\ . To get rid of this purely kinematical dependence, one introduces 
reduced multipoles Aii± via 

= ^f. (20, 

In the isospin-symmetric case, the amplitude for producing a pion with Cartesian isospin 
index a can be decomposed as 

M{ix a ) = x^e a3 VM ( ^ + r a M (0) + 5 a3 M w ) X i, (21) 
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TABLE I: LECs determined from other processes. 



LEC Source 

ci proton mass m p = 938.272 MeV [38] 

C2, C3, C4 pion-nucleon scattering [39] 

C6, C7 magnetic moment of proton (fj, p = 2.793) and neutron (fi n = —1.913) [38 
d±Q axial- vector coupling constant gA = 1.2695 [38] 

c?i8 pion-nucleon coupling constant g n N = 13.21 [41] 



where Xi an d Xf denote the isospinors of the initial and final nucleons, respectively, and r a 
are the Pauli matrices. Here, the only relevant physical channel is given by 

M( 7 + p 7T° +p) = M (0) +M (+) . (22) 

This will be important when it comes to the determination of the LECs in the next section. 



III. THEORETICAL FRAMEWORK 



Chiral perturbation theory is the low-energy effective field theory of QCD based on an 



approximate chiral symmetry and its spontaneous breakdown [3l|, [32 ]_. In the one-nucleon 



sector, pions and nucleons are used as the effective degrees of freedom |33j. Starting from the 



most general Lagrangian in combination with a suitable power-counting scheme, observables 
are calculated in a momentum and quark-mass expansion. At first, the correspondence be- 



tween the chiral and loop expansions known from the mesonic sector seemed to be lost |33 
leading to the development of HBChPT [ll|, 12] in which one projects onto large and small 
components of the nucleon field and, finally, restores a systematic counting scheme. This 
framework was successfully applied to many processes, including pion photo- and electro- 
production [l3 -15, 34 -36]. However, giving up manifest Lorentz covariance may, in certain 
cases, lead to the wrong analytic behavior of loop amplitudes [16J. The seeming power- 
counting problem of RChPT for nucleons was addressed using appropriate renormalization 
schemes such as, e.g., infrared regularization [l6[ or the extended on- mass-shell (EOMS) 
scheme 

The Lagrangian relevant to the one-nucleon sector consists of a purely mesonic part (£„-) 
32] and a part containing the pion-nucleon interaction (£ n N) (H, 3~7l |. 



f _ r(2) , W4) 1 Al) _,_ f\*) 1 /»W , r\^i 



,(2) 



.(3) 



,(4) 



(23) 



The superscripts refer to the chiral order of the Lagrangians and the ellipsis stands for the 
neglected, higher-order contributions. The Lagrangian contains a large number of LECs. 
Their numerical values cannot be derived from the effective field theory itself but are de- 
termined by adjusting them to experimental data. Most of the LECs also enter simpler 
observables such as, e.g., form factors or the pion-nucleon coupling. In Table HI we display 
LECs of £ n N which have been extracted from processes other than pion photoproduction. 
In the following, we focus on the contact interactions specific to pion photoproduction at 
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0(q 3 ) and 0(q A ). The relevant part of the Lagrangian is given by 



^ = 2mV [^ e ^ Tr D P* + H ' c - 

C?9 



+ 



Zffl ]\f 

(4) _ e 48 



[z*e^Tr (/+ + 2wW) + H.c] , 

4^ [^Tr (/+ + 2«W) ^76^* + H 



(24) 



2^ (/J + 2 ^) VfcT^'* + H. 



- £^ [^ A [^a, Tr (/+ + 2^)] 7 57 M ^* + H.c.] 
-^[^Tr(/+^) 757 ^ + H, " 



4m 



-jv L 



MTr (/+ ^) 757"^* + H.c. 



+ 



.c. 



^ - ' "J V 

2^ (/JA*) 757 A ^* + H 

- [zf Tr ( M A [£ A , /+ ]) 757 M ^^ + H.c. 

- i^L [*Tr (/+ + 2«W) X-757 M ^^ + H.c] 



4m 



*Tr(/+X-)757 M ^ + H, 



(25) 



where H.c. refers to the Hermitian conjugate. The nucleon is represented through the 
isospinor field the pion appears after expanding the so-called chiral vielbein «„, and 



the photon is contained in the field-strength tensors /+ 



reader is referred to Ref. 37 



and v$ . 



For further definitions, the 
In neutral pion photoproduction on the proton, only half of 
the LECs listed above can be determined, because the (0) and (+) components of the isospin 
amplitudes [see Eq. (T2TT) ] both contribute in the same way (see Appendix This reduces 
the number of independent LECs from twelve to six. In HBChPT, so far only five constants 
were considered, because that calculation took only S and P waves into account. It turns 
out that at chiral order 0(q A ) another LEC appears which mainly affects the multipole E 2 -. 
In HBChPT one can rearrange the LECs such that two enter Eq + and each of the three P 
waves comes with its own LEC. In the relativistic case the situation is more involved. A 
unique matching of the LECs to the multipoles cannot be done, which can be nicely seen in 
terms of the 1/mjv expansion of the multipoles. At leading order, one reproduces the result 
of Ref. " 
Let 
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At higher order, different linear combinations of the constants appear, 
us now address the renormalization condition. In the EOMS scheme, only terms 
explicitly violating the power counting are subtracted. From the six independent LECs of 
the contact diagrams only one is of 0(q 3 ). The other constants are of 0(g 4 ) and, therefore, 
are not necessary to subtract power-counting-violating contributions, as we calculated the 
process up to and including 0(g). After a heavy-baryon expansion, the combination d$ + dg 
appears only in the multipole P3 at order 0(q 3 ) (see Appendix |A]) and, therefore, power- 
counting- violating terms appear only in this multipole. Furthermore, only diagrams of 0(q A ) 
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can produce such contributions. The standard procedure to access the numerical value of 
such an LEC is through adjustment to experimental data. Here, we exploit this fact to avoid 
the calculation of the power-counting-violating contribution (see Appendix [B] for further 
details). 

Another issue in pion production is isospin symmetry. As we work in the isospin- 
symmetric case, the cusp in the E 0+ multipole cannot appear. In Ref. [36], instead of 
the neutral-pion mass the mass of the 7r + was used in the loops, leading to a phenomenolog- 
ically correct description of the cusp. We also use this method to reproduce the cusp. Even 
though the effect is much smaller, we also use the neutron mass instead of the proton mass 
within the loops. The error one introduces this way is formally of higher order, because the 
mass difference of charged and neutral pions and of proton and neutron is of higher order. 



IV. RESULTS AND DISCUSSION 
A. The RChPT calculation 

Up to and including chiral order four, 20 tree and 85 loop diagrams contribute to the 
more general case of pion electroproduction. The topologies can, e.g., be found in Ref. 
jsj. We calculated all diagrams with the aid of the computer algebra system Mathematica 
and the FeynCalc package |42j. Even with a modern computer program such calculations 
are somewhat cumbersome. Furthermore, using computers always requires control over the 
programs used. In our case crossing symmetry and current conservation provide important 
checks on the results. All our results fulfill these requirements analytically. To evaluate loop 
integrals, we made use of the LoopTools package |43j. 

Using the formulas for the multipole decomposition, we are able to project out any desired 
multipole. For our purposes we calculated S, P, and D waves. The unknown LECs were 
determined via a x 2 fit to the latest MAMI data [21]. These were taken over a much wider 
energy range than ChPT can be applied to. Therefore, we had to determine the best energy 
range for a fit. The problem is that higher-energy data have smaller relative errors, leading 
to an increase of their weight in a fit. We found E l ^ b = 165.8 MeV to be a good maximum 
energy for the fit range. In Fig. [2j we show how the x 2 e d changes if one includes all data 
points up to a certain energy ^ ab . max . For comparison we also provide the reduced x 2 e d °f 
the HBChPT fit. 

It turned out that with the existing data we could only fit five of the six LECs. The 
problematic LEC, e 4 g, strongly influences the E 2 - multipole and, indirectly, also E 0+ . Re- 
sults for these multipoles, taking into account, are shown in Appendix In our final fit, 
we decided to set 649 = and obtained a minimal x 2 e d °f 1-22. 

We estimate the errors of our parameters using the bootstrap method (see, e.g., Ref. 



44j for an introduction). Here, we only briefly outline the idea. If one has a data set 
Y = yi, . . . ,y n of length n, one can create m bootstrap samples Yi, . . . , Y m of length n, 
where m should be a sufficiently large number. These new samples consist of the original 
data points, but randomly chosen. This means that in an arbitrary sample Y^ some data 
points appear twice, three or even more times, while others are neglected. Every sample 
can now be fitted in the same way as the original data. One ends up with m values for 
the parameters. The idea behind the bootstrap is that the standard deviation for each 
parameter is an estimate of its error. Our results for the LECs including an error estimate 
are shown in Table [Til 
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FIG. 2: (color online) x 2 e d as a function of the fitted energy range. The black squares and red dots 
refer to the RBChPT and HBChPT fits, respectively. 



TABLE II: LECs of the contact diagrams. The di are given in units of GeV 2 and the ej in units 
of GeV" 3 . The errors stem from a bootstrap estimate (see text for details). 



LEC Value 





■= d 8 + d 9 


-2.31 ±0.02 


e48 


= e48 + e67 


-3.0 ±0.2 


£49 


.= e 49 ± e 68 





£50 


= e 50 + e 69 


-1.2 ±2.1 


hi 


= e 5 i + e 7 i 


2.3 ±1.1 



em := en2 + eii 3 -4.4 ±2.1 



The graphs for the measured differential cross sections and photon asymmetries are shown 
in Figs. |3] and HI The corresponding multipoles are shown in Figs. |5] and [61 The differential 
cross sections agree nicely with our result in the fitted energy range. For higher energies some 
differences between experiment and our calculation become visible. For the asymmetries the 
overall picture is similar; there, the difference at the highest energies can be traced back to 



P2 and P3 which are a little bit too small compared to the single-energy fit shown in Ref. [21 



The other multipoles agree up to approximately E l ^ h = 170 MeV with the single-energy fits. 
In Fig. |6l the physical P-wave multipoles are shown. Using this representation, one gets a 
clearer picture on the deviations from the data. The multipoles E±+ and Mi_ agree nicely 
with the experiment. For the Mi+ one can see a rising of the data above E\ ah = 170 MeV, 
which is related to the A resonance. As we did not include it explicitly, this rising does not 
appear in our curve. The most important D wave, E 2 -, is discussed in Appendix O 

As we now have the important multipoles in the threshold region, we can make some 
predictions for upcoming experiments. Therefore, in Fig. [7] we show the polarization ob- 
servables T, P, E, E, G, and H. Additionally, we show the predictions of the DMT model 
[22l | . We show the angular distribution at a fixed energy W = 1090 MeV and the energy 
dependence at either 0^ = 90° or = 45°, depending on the approximate extreme value 
of the observables. We find a good agreement between RChPT and DMT. 
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FIG. 3: (color online) Differential cross sections <7o(©7r) in /ib/sr as a function of the cm production 
angle Q n . The graphs are shown for increasing pion energies in the cm frame/photon energies in 
the lab frame, both given in units of MeV. The solid (black) curves show the results in RChPT at 
0(q 4 ), the dashed (red) curves show the same chiral order in HBChPT. The fits make use of data 



up to and including E, 



lab 



165.8 MeV, i.e., the first nine figures. The data are taken from Ref. 
11 




FIG. 4: (color online) Photon asymmetries £ as a function of the cm production angle 0^. The 
graphs are shown for increasing pion energies in the cm frame/photon energies in the lab frame, 
both given in units of MeV. The solid (black) curves show the results in RChPT at 0(q ), the 
dashed (red) curves show the same chiral order in HBChPT. The fits make use of data up to and 
including Elf 3 = 165.8 MeV, i.e., the first nine figures. The data are taken from Ref. [2ll |. 
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FIG. 5: (color online) Real parts of the S- and P-wave multipoles as a function of . The solid 
(black) curves show the result in RChPT at 0(g 4 ), the dash-dotted (red) curves show the same 
chiral order in HBChPT. The dashed (black) lines show the RChPT result without the LECs, the 
dotted (red) lines show the HBChPT result without the LECs. The long-dashed (green) curves 
show the HBChPT fit of Ref. 



211 ] . The data points are taken from a single-energy fit from Ref. 



B. Comparison with HBChPT 



The HBChPT calculation of Ref. 
new experimental data 14, 15, 29, [30|, |45|, 146 



151 has been analyzed several times in the light of 
Here, we fit the LECs in the same way to 
the experiment as in the relativistic case, allowing for a better comparison between both 
calculations. Note that the HBChPT fit of Ref. H extends to a larger value of E^ h 
which, as we have checked within our calculation, partly explains the difference between the 
two HBChPT results in Fig. [5] Another source for the discrepancy is the use of different 
values for the coupling constants and different fitting procedures. 

As mentioned before, D waves are expected to be of some importance for extracting the 
E 0+ multipole. We do not have the heavy-baryon result for the D waves. Therefore, we use 
the Born terms to calculate D waves. In addition, we included the heavy-baryon-expanded 
contribution of the sixth LEC to E2-. Eventually, we experienced similar problems fitting 
this LEC as in the RChPT case. For that reason, we also set this constant to zero and used 
only the Born terms for the D waves. With a value of 1.11, the Xred * s better than in the 
relativistic case. One can see this in the observables too, as the heavy-baryon calculation 
seems to agree slightly better with the data. The multipoles also support this picture. The 
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FIG. 6: (color online) Physical P-wave multipoles as a function of . The solid (black) curves 
show the results in RChPT at 0(q 4 ), the dashed (red) curves show the same chiral order in 
HBChPT. The data points stem from a single-energy fit from Ref. 211 ] . 



TABLE III: Values of the LECs in HBChPT obtained from a fit to Ref. Q. We decided to set 
e&9HB = (see text). The errors stem from a bootstrap estimate (see text for details). 

~~ o'i (15.2 ± 2.7) GeV" 4 

a 2 (-7.6 ± 2.5) GeV" 4 
£i 33.3 ±0.5 

£ 2 -31.8 ±0.5 

b p (20.9 ± 0.1) GeV" 3 

e A9HB 



LECs we obtained are listed in Table IHIl On the other hand, the convergence properties of 
the relativistic result look more favorable. To illustrate this we also display the multipoles 
one gets, when switching off the LECs. The difference between the final result and this 
case is some indication of how good the series converges. Comparing the size of the LECs 
leads to the same conclusion. The most dramatic effect appears in the case of P3, where in 
HBChPT the LEC completely dominates the Born and loop terms. This gives us confidence 
that certain higher-order terms are very important here. The relativistic calculation keeps 
some corrections up to infinite order and this improves the convergence behavior. 
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FIG. 7: (color online) The polarization observables T, P, E, F, G, and H. The left column shows 
the angular dependence at W = 1090 MeV and the right column shows the energy dependence of 
T, E, and G for Q n = 90° and of P, F, and H for Q n = 45°, respectively. The solid (black) curves 
show the RChPT results and the dashed (blue) curves show the DMT model. 
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Another interesting quantity one can derive from pion photoproduction is the so-called /3 
parameter 30] of the unitarity cusp 23[. It is linked to pion-nucleon scattering and charged 
pion photoproduction via 



/3 = M, + Re[Eo + (7,vr + )]a c 



7r n 



7c°p) = (3.43 ±0.08) x 10~ 3 /M n+ , 



(26) 



where the numerical estimate [2l| makes use of isospin symmetry to replace a ccx (7r + ri — > n°p) 
in terms of the experimentally known scattering length a cex (iT~p — > 7r°n). Close to threshold, 
unitarity connects this parameter to the imaginary part of E 0+ (j, 7T°), 

(27) 



(28) 



Im [E 0+ (7,7r )] =(3q + , 
where q + is proportional to the three- momentum \q + \ of a ir + in the cm frame, 

q + = \q+\/M^+. 



To pin down the numerical value of j3, we fit the imaginary part of Eq + to the following 
series: 

^lab ^lab,thr 



Im [E 0+ (j, 



7T 



(29) 



In case of the relativistic calculation, we get /3r = 3.16 x 10 3 /M 7r + and 7^ = —1.08 x 
10- 3 / M ^+> and HBChPT results in f3 HB = 2.83 x IQ~ 3 /M n + and ^ HB = -1.97 x IQ~ 3 /M n +. 
Both results are predictions, as all LECs were fixed in other processes, including pion-nucleon 
scattering. Nevertheless, both results are too small compared to the experimental value of 
(3.43 ± 0.08) x 10~ 3 /M 7r + [2lJ. The relativistic result is somewhat closer, indicating that, 
again, certain higher-order contributions are important. In Fig. [HI the imaginary part of 
Eq + is shown function of E 1 ^ . 
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FIG. 8: (color online) Imaginary part of £^0+ as a function of . The solid (black) curve shows 
the result in RChPT at 0(q 4 ), the dashed (red) curve shows the same chiral order in HBChPT. 
The (green) band shows the result from unitarity with /3 = (3.43 ± 0.08) x 10 _3 /M 7r + and 7 = 0. 
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V. SUMMARY 



In this article we presented a full RChPT calculation up to 0(g 4 ) and one-loop order for 7r° 
photoproduction on the proton. The amplitude was calculated with the aid of Mathematica 
and additional packages. Several tests were utilized to check our results. 

The next step was an analysis of the results in terms of multipoles. We took S, P, and 
D waves into account, as these are the only relevant multipoles in the threshold region. 
The LECs of the contact diagrams were fitted using data of the latest MAMI experiment 
for differential cross sections and photon asymmetries. Our calculation agrees well with the 
experimental data up to photon energies of approximately 170 MeV in the lab frame. 

We also discussed some of the properties of the LECs concerning the multipoles. We 
found that there are six LECs, of which one (649) is very important for the determination 
of E 2 - and, therefore, also for E 0+ . With the existing data we cannot narrow down this 
constant, because an unconstrained fit gives an unnaturally large value for this LEC. Hence, 
we decided to neglect it in our main analysis. With the multipoles at hand we also gave 
some predictions for polarization observables for future experiments. 

In addition, we compared our results with HBChPT. Even though the latter seems to 
describe the experimental data slightly better, the corresponding LECs have rather large 
values. This indicates that the relativistic calculation contains important higher-order ef- 
fects. This can also be seen in the /3 parameter, where the RChPT result is closer to the 
commonly accepted value stemming from unitarity. 



Our numerical results are available through a web interface |47J]. The freedom to change 



the LECs gives one the opportunity to explore our results further in the light of new upcom- 
ing data. We also give results for pion electroproduction and for all four physical reaction 
channels. The details will be discussed in a forthcoming publication |48[. 
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Appendix A: Heavy-baryon expansion of the contact contributions 

Here, we provide the results of expanding the contact contributions to Eq + , the three P 
waves, and E2- in powers of l/rajy up to and including next-to- leading order (heavy-baryon 
expansion) : 

M = M m + — M [1] + I 4rV ( A1 ) 
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Using the definitions of Table [Hi we obtain at leading order: 



p [0] _ e(6e 48 + 2e 49 - 4e 50 + 3e 5 i)E* e (3e U2 + e 49 ) M 2 E 7T 

ho+ ~ UtTF 6t7f ' {AZ) 

p[o] _ e(2g 48 + hi)El 

P n = ^ (A4) 

if = (A5) 

At next-to-leading order, the results read 

[i] e(3en 2 + e 49 )M 4 e(6e m - 12e 48 - 4g 49 + 12g 50 - 5e 5 i)M 2 £ 2 
0+ 12vrF 24ttF 

e(-3e 48 + e 49 - 2e 50 - 2e B i)S* , , 

+ 127F ' (A7) 
p[i] e(2g 48 - 2g 50 + e 5i )^ e(-g 112 + 2g 48 + e 51 )M 2 E w 

1 ~ AttF AttF ' 1 ' 
= [!] e(6g 49 + 4g 50 -3g 5 i)£ 3 e(e m - 2g 48 )M 2 ^ 

P * " 24^F + 4^F ' (A9) 

if = ^ (A10) 

3 2ttF ' v ; 

= [i] _ ee 49 M 2 _ e(6g 48 + 4g 49 + g 5 i)E 2 

2 ~ 12vrF 48vrF ' [ ' 

As one can clearly see, only the lowest order allows for a rearrangement of the LECs such 
that one can uniquely assign them to the different multipoles. For P 3 this also works up to 
and including first order. However, the other multipoles generate new mixings of the LECs 
at next-to-leading order. 



Appendix B: Renormalization of power-counting-violating contributions 

To some extent the EOMS scheme can be utilized to renormalize diagrams without ex- 
plicitly calculating the power-counting-violating part of the diagrams. Here, we explain this 
statement using a generic example, namely, the mass of a particle. Let us assume for the sake 
of simplicity that, after renormalization, power counting predicts a tree-level contribution 
of chiral order 0(q°) and a loop contribution of chiral order 0(q 2 ). Before renormalization, 
the mass is of the form 

m = LEC° + Loop , (Bl) 

where LEC° represents an unknown bare LEC and Loop represents the unrenormalized 
loop contribution. In the following, we neglect any ultraviolet divergences in these expres- 
sions, i.e., we assume that they have been taken care of by applying the modified minimal 
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subtraction scheme of ChPT [33|. We indicate this fact in terms of a superscript r. The 
infrared regular part of the loop contribution (Loop r IR ) can be symbolically written as 

Loopj R = «q + a r 2 q 2 + a\q A H , (B2) 

where q is a small quantity. For notational convenience, we have assumed only even powers 
of q. According to the above assumption, a r Q violates the power counting. Renormalizing 
the loop contribution in the EOMS scheme amounts to subtracting the power-counting- 



violating term a r from Loop^ [17] . In other words, Loopj R is replaced by Loopf^, and Eq. 
(IB 2 1) simply becomes 

Loopf fl = a R q 2 + afq 4 + ■■■ . (B3) 

Note that, in general, the higher-order coefficients 012, «4, • • • are expressed in terms of EOMS- 
renormalized quantities. In this case, the expression for the mass reads 

m = LEC R + Loop 75 , (B4) 

where the renormalized constant LEC R now absorbs the power-counting- violating part. 
Comparing Eqs. ( 1B1[) and ( 1B4I) . the following connection can be made 

LEC r + a r = LEC R . (B5) 

The two expressions LEC r and LEC R only differ by the (numerical) value of a r . In con- 
clusion, adjusting LECs numerically is sufficient to renormalize diagrams. However, this 
procedure also has some drawbacks. First of all, only the sum of the adjusted LECs and 
the loop diagrams satisfies the power counting whereas in the standard EOMS procedure 
each renormalized diagram satisfies the power counting separately. As a consequence, in 
our calculation we cannot separate the fourth-order loop correction for P 3 . Another general 
problem is that the LEC does not necessarily have to be of natural size anymore, as it 
contains a power-counting- violating part. In our specific case this is true for loop diagrams 
at fourth order that are renormalized in terms of dg. Nevertheless, in the present case, it 
turns out that this part is either small or has the opposite sign and same magnitude of 
the numerical contribution of the renormalized coupling, because even though dg contains a 
power-counting-violating part it is of natural size. 



Appendix C: The LEC 649 

In principle, the LEC 649 appears in all multipoles, but as one can see from Eq. (IA6j) it 
mainly affects E2-. This multipole mixes strongly with Eq + . When performing a completely 
unconstrained fit, the solution with the lowest x 2 e d now * s Xred = 1-14 and changes the values 
of E2- significantly. The numerical value for e^g becomes unnaturally large, as one can see 
from Table IIVI The multipoles affected most are shown in Fig. EO From the point of view 
of naturalness, it is very unlikely that an LEC picks up such a large value. Therefore, we 
believe such a large value to be an artifact until beam-target double polarization observables 
E and F are available. This artifact can also be seen in the heavy-baryon calculation. There, 
we used the lowest-order contribution of the LEC for E 2 - [see Eq. (1A6j) ]. as it would appear 
in exactly the same way in a true heavy-baryon calculation at 0(q A ) including D waves. 
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TABLE IV: LECs of the contact diagrams in a fit including 649. The di are given in units of GeV 
and the in units of GeV -3 . 
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FIG. 9: Results for the E$ + and E2- multipoles with 649 = (solid line) and 649 = —35.7 GeV 
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